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(57) Abstract: The invention relates to processing seismic 
data that contains both a desired signal and swell noise. The 
method is applicable to seismic data is in the frequency-space 
domain method and comprises determining a signal-only pre- 
diction filter from the seismic data at a first frequency at which 
swell noise is not present, and applying the prediction filter 
to seismic data at a second frequency at which swell noise 
is present. This attenuates swell noise in the seismic data at 
the second frequency. In one embodiment, the prediction fil- 
ter for a frequency f L , the lowest frequency at which swell 
noise is present, is determined from seismic data at frequency 
fH+i» where f H is the highest frequency at which swell noise 
is present. The prediction filter for a frequency f^ is deter- 
mined from seismic data at frequency f H+ 2» and so on, so that 
the prediction filter for frequency f H+ i is determined from seis- 
mic data at frequency f2H-L+i* 1° tn * s wav » tne prediction fil- 
ter for each frequency in the swell noise range is determined 
from seismic data outside the swell noise range, but at a fre- 
quency close to the swell noise range, and this provides effec- 
tive extraction of the signal and attenuation of die swell noise. 
Furthermore, no prediction filters is applied to data of more 
than one frequency, so that an error in the determination of a 
particular prediction filter will have only limited effect on the 
processed data. 
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Processing Seismic Data 



The present invention relates to processing seismic data and, in particular, to processing 
marine seismic data that includes swell noise so as to attenuate or remove the swell 
noise. 

In a seismic survey, acoustic energy is emitted from a source and is detected by a sensor 
(generally referred to as a "receiver") located at a distance from the source. Some of the 
acoustic energy emitted by the source passes into the earth's interior and is reflected by 
geological structures within the earth. Information about the geological structure of the 
earth's interior can be derived from the reflected acoustic energy incident on the 
receiver. 

In a marine seismic survey both the acoustic source and receiver are disposed within a 
water column, such as the sea. Figure 1 is a schematic illustration of one marine 
seismic survey in which an acoustic energy source 1, such as an airgun, is suspended in 
the water column 2 from a survey vessel 3. An array of receivers 5 such as, for 
example, hydrophones, is also disposed in the water column 2, and the receivers 5 are 
suspended from a float 6 that is towed by a second survey vessel 4. Such a receiver 
array is known generally as a 'towed marine" receiver array. 

When the acoustic source 1 is actuated it emits acoustic energy. Acoustic energy 
emitted by the source propagates generally downwards through the water column, 
passes into the earth's interior 7 and continues to propagate generally downwards 
through the earth's interior until it undergoes partial reflection at a geological structure 
within the earth's interior. Reflected acoustic energy propagates generally upwards 
through the earth's interior 7 and through the water column 2 and is incident on the 
array of receivers 5. Seismic data acquired by the receivers 5 following actuation of the 
acoustic source 1 are transmitted to the second survey vessel 4, for example by an 
electrical or optical link. The data may be analysed in a processor 8 located on the 
second survey vessel or may be stored for subsequent analysis. The principle of seismic 



WO 2004/051315 



PCT/GB2003/005278 



surveying is to obtain information about te earth's interior by suitable processing of the 
seismic data acquired at the receivers 5. 

In general, more than one geological structure within the earth's interior acts as a partial 
reflector for acoustic energy. The seismic data acquired at a receiver 5 foUowmg 
actuation of the acoustic source 1 therefore contain a number of spatially coherent 
events, each of which corresponds to reflection of acoustic energy at an associated 
geological feature wimin the earth that acts as a partial reflector of acoustic energy. 

One problem in processing towed marine seismic data is that data acquired in rough sea 
conditions may contain swell noise. Swell noise may be either coherent or random, and 
occurs at low frequencies, generally below 30Hz. Swell noise generated by a rough sea 
interferes destructively with the desired events in the seismic data (i.e„ the events 
arising from reflection of seismic energy by geological features within the earth). Swell 
noise must be removed from the seismic data at some stage during the processing of the 
data, and the removal of swell noise is one of me most important steps leading to final 
interpretation of the acquired data and subsequent mapping of the earth's interior. 

A number of prior art methods for attenuation of swell noise are known. 

GA McMechan has proposed, in "Automatic editing of noisy seismic data", 
Geophysical Prospecting, Vol. 37, PP 875-892 (1989), a method in which me low 
frequency energy content of a seismic trace is evaluated. A trace is discarded if its low 
frequency energy content exceeds a given threshold. This method has the severe 
disadvantage that the desired signal in me trace is also lost in me process. 

H Truman Holcombe and Robert S. Wojslaw have proposed, in "Spatially weighted 
trim stacking, a technique for prestack suppression", Society of Exploration 
Geophysicists, PP 1157-1160 (1992), a method that involves locating swell noise m a 
trace and replacing the affected portion of the trace with an interpolation which makes 
use of the neighbouring traces. Since mis method requires an interpolation usmg 
neighbouring traces reflector dip information is merefore required, so that this method is 
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a costly and often inaccurate process. Furthermore the dip search can only provide the 
dominant dip. 

Guillaume Cambois and Jacques Frelet have proposed, in "Can we surgically remove 
swell noise ?", Society of Exploration Geophysicists, 1995, a noise free prediction filter 
that is designed and applied in the f-x (frequency-space) domain. This filter is quite 
efficient at filtering out swell noise when it is random, but coherent swell noise is not 
eliminated by this melhod but rather is enhanced. This "one step ahead" prediction 
filter is a coherent enhancement process and hence cannot separate coherent swell noise 
from a coherent signal 

Ali Ozbek has proposed, in a patent application entitled "Adaptive Seismic Noise and 
Interference Attenuation Method" an adaptive beaniforming approach for attenuating 
the various types of coherent noise encountered in seismic data acquisition. This method 
applies an adaptive (t-x) time space filter in those parts of the frequency-wavenumber 
space that contain the noise. Seismic signals are passed while coherent noise and 
interference components that are temporally and spatially not stationary are adaptively 
filtered This method can be quite effective but can be computationally intensive. 

til 

L.Canales has also proposed a method of attenuating swell noise, in the SGE 54 
Annual International Meeting (1984). This is a "one step ahead" wiener prediction 
filtering process carried out in the f-x (frequency space) domain. Coherent signal and 
coherent swell noise are both enhanced. 

S.Spitz has proposed, in "Adaptive Dip Filtering in the F-X Domain", SEG 63 rf Annual 
International Meeting (1993) a method in which the prediction filters corresponding to 
all the dips within the data are evaluated through polynomial factorisation of the fall 
prediction filter. All the dipping events, coherent noise include, within Ihe input data are 
men solved for using a least square approach. This is a difficult task as all the dips need 
to be accounted for. 
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Proprietary modules that provide similar functions are known, and these generally 
involve comparing a time window within a given trace with surrounding traces to check 
for the presence of noise, usually against a user-specified threshold value. Noise above 
the specified threshold is removed. Such methods are, however, very sensitive to the 
threshold, and are likely either to remove too much from the data including some of the 
signal of interest or to leave some residual noise. 

The present invention provides a method of processing seismic data, the data being in 
the frequency-space domain, the method comprising the steps of; determining a 
prediction filter from the seismic data at a first frequency at which swell noise is not 
present; and applying the prediction filter to seismic data at a second frequency at which 
swell noise is present, thereby to attenuate swell noise in the seismic data at the second 
frequency. 

The prediction filter separates the swell noise from the signal of interest in the data, and 
attenuates the swell noise while fully preserving the signal of interest The invention 
provides effective attenuation of both coherent and random swell noise without 
affecting the signal of interest. 

A second aspect of the invention provides a method of marine seismic surveying 
comprising actuating a source to emit acoustic energy; acquiring seismic data in the 
frequency domain; and processing the data according to a method as defined above. 

A third aspect of the invention provides an apparatus for processing seismic data, the 
data being in the frequency-space domain, the apparatus comprising means for 
determining a prediction filter from the seismic data at a first frequency at which swell 
noise is not present; and means for applying the prediction filter to seismic data at a 
second frequency at which swell noise is present, thereby to attenuate swell noise in the 
seismic data at the second frequency. 

The apparatus may contain a programmable data processor. 
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A fourth aspect of the invention seismic surveying arrangement comprising: a source of 
acoustic energy; means for acquiring seismic data in the frequency-space domain; and 
an apparatus as defined above. 

A fifth aspect of the invention provides a storage medium comprising a program for a 
data processor of an apparatus as defined in above, and a sixth aspect of the invention 
provides a storage medium containing a program for controlling a data processor to 
perform a method as defined above. 

Preferred features of the invention are defined in the remaining claims. 

A preferred embodiment of the invention will now be described in detail, by way of an 
illustrative example, with reference to the accompanying figures in which; 

Figure 1 is a schematic illustration of a towed marine seismic survey; 

Figure 2 illustrates a typical seismic data trace, in the frequency domain, obtained by the 
seismic survey of figure 1; 

Figure 3 is a flow diagram of a method according to one embodiment of the present 
invention; 

Figures 4(a) to 4(c) illustrates the results of the invention on synthetic seismic data; 

Figures 5(a) to 5(c) illustrates the results of the invention as applied to single sensor 
seismic data; 

Figure 6 is a schematic block diagram of an apparatus according to the present 
invention. 

A seismic data trace acquired at one of the receivers 5 in the seismic surveying 
arrangement of figure 1 records the amplitude of a seismic parameter, such as the 
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pressure or a component of particle motion, as a function of time. Examples of typical 
seismic data traces are shown in Figure 4(a) and 5(a) (these traces include swell noise) 
and in Figures 4(b) and 5(b) (which show the traces of Figures 4(a) and 5(a) after 
removal of swell noise). The acoustic source is actuated at time t = 0, so mat a trace 
measures the variation in amplitude of the parameter following the actuation of the 
acoustic source. It will be seen from Figures 4(b) and 5(b) that, when the acoustic 
source is actuated, there is initially, in the absence of noise, no signal acquired at the 
receiver. The first event in the seismic data trace generally arises from acoustic energy 
that has travelled direct from the source 1 to the receiver 5 through the water column, 
without entering the earth's interior. The time delay between the actuation of the 
acoustic source and this event (which is know as the "direct evenf ') depends on the 
speed of propagation of acoustic energy through the water column and the distance 
between the acoustic source and the receiver. The direct event is then followed by 
events that arise from acoustic energy that has entered the earth's interior 7 and has 
undergone reflection at a geological feature within the earth. 

The seismic data trace of figures 4(a) to 5(c) are shown as continuous lines. In practice, 
however, a receiver will generally sample the seismic parameter with a sampling 
frequency f s (that is, the parameter is sampled at time intervals of l/f s ). The sampling 
interval l/f s may be, for example, 1 , 2 or 4ms. A data trace therefore consists of a series 
of discrete values for the parameter, one obtained hi each sampling operation. 

The seismic data trace in the time domain may be Fourier transformed to produce a 
record of the amplitude of Ihe seismic parameter as a function of frequency. Figure 2 is 
a schematic illustration of a typical data trace in the frequency domain, as would be 
obtained by transforming one of the data traces of Figure 4(a) or 5(a) to the frequency- 
space domain by, for example, using a Fourier transform. Since the data trace in the 
time domain is obtained by a sampling process, with a sampling frequency fs, the 
transformed data trace of figure 2 extends only to a maximum frequency of fey* where 
f^ is related to the sampling frequency fs by the Nyquist theorem The data trace is 
shown as a continuous trace in Figure 2 for convenience. As noted above, however, a 
seismic data trace in the time domain in practice consists of a series of discrete values 
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rather than a continuous trace. The data trace in the frequency-space domain of figure 2 
will therefore also consist of a series of values at discrete frequencies fj, where the index 
j extends to the highest value J at which fj does not exceed the Nyquist frequency. 

Swell noise is a low-frequency noise, typically having a maximum frequency of around 
30 Hz. The frequency-space domain data trace of figure 2 therefore contains two 
regions - a low frequency region 9 that contains both the signal of interest and swell 
noise, and a higher frequency region 10 that contains only the signal of interest In 
figure 2 the frequency region 9 that contains both swell noise and the signal of interest 
extends from frequency f L to frequency fe, and that the frequency band that contains 
only the signal of interest is in the range from frequency f H +i to f^q. 

The lower frequency limit fL, of the band containing swell noise will generally be zero. 
If f L is not zero, the data at frequencies below ft. represents a combination of signal and 
noise other than swell noise. 

In the method of the invention, a signal prediction filter is determined from the seismic 
data trace of figure 2 at a first frequency at which swell noise is not present - that is, at a 
frequency of fe+i or greater. The prediction filter is then applied to seismic data at a 
second frequency at which swell noise is present - that is, at a frequency between f L and 
f H - so as to attenuate swell noise in the seismic data at the second frequency. The 
method is preferably repeated to obtain prediction filters for one or more other 
frequencies at which swell noise is present in the seismic data, for example by deriving 
a filter from data at a third frequency (at which swell noise is not present) and using this 
prediction filter to attenuate swell noise in the seismic data at a fourth, frequency (at 
which swell noise is present). In a particularly preferred embodiment, an associated 
prediction filter is determined for each frequency in the swell noise range from fi. to f H , 
to enable swell noise to be attenuated from every frequency value in the swell noise 
range. 



In the ideal case, the prediction filter that was to be applied at a frequency fj would be 
derived from seismic data for that frequency. In the present case, however, this ideal 
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criterion would require that the prediction filter to be applied to seismic data at 
frequency in the range from f L to f H would be derived from the seismic data for that 
frequency — but this would mean deriving the prediction filter from seismic data that 
contains swell noise in addition to the signal of interest According to the invention, 
therefore, the prediction filter for a frequency component fj is derived from the seismic 
data for a frequency that is close to j§ but that is outside the frequency range within 
which swell noise occurs. This is a reasonable approximation, provided that the 
characteristics of the signal of interest are identical, or very similar, for the frequency fj 
and for the frequency at which the filter is derived. 

In principle, a single prediction filter determined from data for a single frequency at 
which swell noise is not present could be applied at every frequency at which swell 
noise is present This however has the potential disadvantage that if the filter is 
corrupted for any reason, then the output obtained by applying the filter will be noisy at 
every frequency in the swell noise range. It is therefore preferable that a plurality of 
prediction filters, each determined from data at a different frequency, are used so that 
the effect of one corrupted filter will be reduced. It is particularly preferable that each 
prediction filter is applied to data at only one frequency in the swell noise range. This 
means that the number of different prediction filters required is equal to the number of 
frequencies at which swell noise is to be attenuated - and, as noted above, each filter is 
preferably determined at a frequency close to the swell noise range. 

In a particularly preferred embodiment of the invention, the filter to be applied to the 
data at frequency ft is derived from data at frequency f H +i. The frequency f H +i is, as 
noted above, outside the range at which swell noise occurs, and is the lowest frequency 
outside the range at which swell noise occurs. The signal characteristics at frequency 
fk+i should be similar to those for frequency f L . 

In this embodiment, the filter to be applied at frequency ft+i is derived from the seismic 
data at f H +2. Again, the frequency f H +2 is outside the frequency band within which swell 
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noise occurs, and the signal of interest at frequency f H+2 should be similar to the signal 
of interest at fL+i. 

Furthermore, in this embodiment, for each frequency within the frequency band for 
which swell noise occurs, the associated prediction filter is determined from the seismic 
data at a frequency fk«m- Thus, finally, the prediction filter for the highest frequency 
f H in the swell noise band is derived from me seismic data at frequency fai-L+i- 

The frequency range in which the prediction filter is derived thus starts at fe+i, and this 
frequency fulfils two important criteria. Firstly, data acquired at this frequency, or at a 
higher frequency, is assumed to be free of swell noise. Secondly, mis frequency is as 
close as possible to frequency f L , so that the signal of interest should be similar at 
frequencies fen and f L . Furthermore, in this embodiment each filter is applied only to 
data at one frequency, so that the effect of a corrupted filter is nmiimised. 

Figure 3 is a block flow diagram of one embodiment of the present invention. 

Initially at step 1 seismic data are obtained, for example, using a towed marine seismic 
surveying arrangement of the type shown in figure 1. The seismic data will consist of a 
number of traces generally similar to those shown in figure 4(a) or 5(a). 

The invention may alternatively be applied to pre-existing seismic data. In this case, 
step 1 is replaced by step 2 of retrieving pre-existing seismic data from storage. 

At step 3, the seismic data are Fourier transformed to produce seismic data traces of the 
general type shown in figure 2 that contain information as to the frequency dependence 
of the amplitude of a seismic parameter. 



It should be noted that, in an embodiment where stored seismic data are retrieved from 
storage in step 2, step 3 may be omitted if the retrieved seismic data are already in the 
frequency domain; this is indicated by a broken line. 



2004/051315 



PCT/GB2003/005278 



10 



At step 4, a frequency counter i is initialised, and at step 5 a filter F s is determined for 
an initial frequency component fi in the range where swell noise is present in the seismic 
data. Conveniently, step 5 is initially performed for me lowest frequency component at 
which swell noise is present, that is, for frequency f, = f L (and in this case step 4 will 
have set the counter i equal to L). As mentioned above, in one preferred embodiment 
me prediction filter for frequency f L is determined from seismic data at frequency f H +i. 

At step 6, the data at frequency for all traces, are operated on by Ihe filter F s 
determined for frequency ft in step 5, according to equation 1. 

At step 7, it is determined whether the filter has been applied a pre-determined number 
(n) times. The pre-determined number n may be specified by the user for processing a 
particular data set If step 7 yields a "no" determination, at step 8 the filter F s is re- 
applied to the results obtained in step 6. Step 7 is then repeated, and if a "yes" 
determination is still not obtained, steps 8 and 7 are repeated until a "yes" determination 
is obtained 

For example, a user could specify that the filter is to be applied three times (n = 3), in 
which case step 8 would be carried out twice (so mat the filter was applied once to the 
initial data and twice to the output of the previous iterations.) If desired, the method 
may comprise the step (not shown) of a user defining the number of times the filter 
should be applied. 

In an alternative embodiment, step 7 is replaced by a quality control step (not shown), 
so that the filter is applied until a desired result is obtained, and is not applied a pre- 
determined number of times. For example, the residual obtained in one step may be 
compared with the previous residual. 

At step 9, it is determined whether the counter indicates that the process has been 
carried out for all frequencies in the range where swell noise is present If the value of 
the counter i is not equal to H so f, is not equal to fa, step 9 yields a "no" determination, 
the frequency counter i is incremented at step 10, and steps 5-8 are then repeated for the 
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next frequency - in this example, f L+ i- Assuming mat a "no" determination is still 
obtained at step 9, the counter is again incremented at step 10, and steps 5-9 are 
repeated until a "yes" determination is obtained in step 9. A "yes" determination will be 
obtained when a method has been carried out in all frequencies in me range ft, to f H - 
that is, for all frequencies at which the data contains swell noise. 

Alternatively, the method may comprises performing steps 5 to 8 initially for fk, in 
which case step 10 comprises decrementing the counter. In this variation a "yes" 
determination is obtained at step 9 once the frequency has been reduced to ft. 

Once a "yes" determination is obtained at step 9, the results of the method may be 
output if desired (not shown). The results output may include, for each original data 
trace, a corrected version of the data trace from which swell noise has been removed. 
The corrected data trace may be obtained, for example, by subtracting the final residual 
obtained at the end of the iteration step process 7 and 8 from the original data trace. 
Additionally or alternatively, the results output may include the residual obtained from 
each data trace. 

Additionally or alternatively the corrected data traces may be subjected to further 
processing steps, to obtain information about the geological structure of the earth's 
interior, and mis is indicated at step 12. If necessary, the corrected data traces are 
transformed back to the time domain at step 1 1 before step 12 is applied These further 
data processing steps are conventional, and will not be described here further. 

Figures 4(a) to 4(c) and 5(a) to 5(c) illustrate results of the present invention. Figures 
4(a) to 4(c) illustrate the invention as applied to synthetic seismic data, and figures 5(a) 
to 5(c) illustrate the invention as applied to actual marine seismic data. 

Figure 4(a) shows synthetic marine seismic data obtained using a simulation method 
that generates swell noise. The data were simulated for a seismic surveying 
arrangement in which the direct arrival occurs approximately 0.8 seconds after the 
acoustic source has been actuated. Figure 4(a) shows approximately 75 simulated data 
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traces, in the time domain. Each trace shows the amplitude (horizontal axis) of the 
pressure as a function of the time (vertical axis) after the actuation of the acoustic 
source. It will be noted that the data traces contain significant variations in amplitude at 
times before the direct arrival. 



Figure 4(b) shows the traces of figure 4(a) after swell noise has been attenuated using a 
method of the invention, and figure 4(c) shows the residual which represents swell 
noise. (It should be noted that although the method of the invention is carried out in 1he 
frequency domain, the results for the corrected data Iraces and the swell noise have been 
transformed back to the time domain to enable easy comparison with the original data 
traces.) 

If figure 4(b) is examined, it will be seen that there is very little structure in the seismic 
data traces before the direct arrival occurs at approximately 0.8 seconds after the 
actuation of the acoustic source. This indicates that the invention is effective in 
removing swell noise since, in the absence of noise, the receivers should not detect any 
acoustic energy from the source before the direct arrival. It will also be noted that there 
is significantly less structure in the seismic data traces for times greater man 
approximately 3 seconds after the actuation of the acoustic source. This indicates that 
most of the structure seen in the original data traces of figure 4(a) at 3 seconds or 
greater after the actuation of the acoustic source does not result from reflection by a 
geological structure but arises from swell noise. 

It will further be noted that the events in the seismic data traces occurring at times from 
approximately 1.6 to 3 seconds after the actuation of the seismic source are much 
clearer in the corrected seismic data traces of figure 4(b) than in the original data traces 
of figure 4(a). This indicates that the attenuation of swell noise has been effective 
throughout the entire data trace. As a result, information about the structure of the 
earth's interior derived from analysis of these events should be more accurate than if the 
original uncorrected data traces of figure 4(a) had been used. 
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Figures 5(a) to 5(c) again show results of the present invention, but applied to actual 
seismic data. These seismic data we acquired using a marine seismic surveying 
arrangement in which the direct arrival occurs approximately 1.8 seconds after the 
actuation of the acoustic source. The figures show data for approximately 70 receivers, 
and these are arranged as a linear array extending generally away from the source. As a 
result, the arrival time of the direct event increases gradually from one trace to the next, 
as a consequence of the increasing distance between the source and the receiver used to 
acquire a particular trace. 

As is clearly seen in figure 5(a), the acquired seismic data contain significant amplitude 
variations for times before the arrival of the direct event and, as explained above, these 
amplitude variations arise primarily from swell noise. 

Figure 5(b) shows the traces of figures 5(a) after attenuation of swell noise using a 
method of the present invention. It will be seen that the amplitude variations in the 
traces for times before the arrival of the direct event have been significantly reduced. It 
will also be seen that many of the events in the traces for long arrival times have also 
been eliminated, and this indicates that these events arise from swell noise rather from 
reflection by a geological feature within the earth. The events of interest in the data 
traces, which generally occur at times between 1.8 and approximately 5-6 seconds after 
the actuation of the source are much clearer in figure 5(b) than in the original traces of 
5(a). 

Figure 5(c) again shows the residual, which represents swell noise. The results for the 
corrected data traces and the swell noise have again been transformed back to the time 
domain to enable easy comparison with the original data traces. 

A preferred embodiment of the filtering method will now be described. 

Let D, denote the input data. This is assumed to consist of N traces identified by the 
index i where 1 <; i <. N. These traces are acquired at the same time but by different 
receivers in a receiver array. Each trace will contain a signal of interest, and will also 
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contain swell noise. The signal of interest in the trace will be denoted by X si and the 
swell noise value in the I th trace will be denoted by X m . Thus: 

D 3 =X S3 +X m 



The signal prediction filter is denoted by F s = [F SI F S2 F^ . . . FgJ, where Fs* 
denotes the A* filter component F s is a general term for a prediction filter derived at a 
particular frequency outside the swell noise range. The index k satisfies 0 £ k < P, 
where P is the number of signal related dips. The filter length is preferably between 3 
and 10 traces. The prediction filter may be derived by any suitable technique such as, 
for example, Wiener filtering. As noted above, in a preferred embodiment the 
prediction filter for frequency ft, is determined from data at frequency fk+i, the 
prediction filter for frequency f L +i is determined from data at frequency f H +2, and so on. 

An appropriate filter F s is applied to each seismic data trace at each frequency in the 
range (ft, - fti) for all traces. The resultant signal is predictable as a function of space 
(Canales 1984, above). 

Filtering the input data with F s only will result in the following system of linear 
equations: 
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The right hand side of equation (1) contains two terms. The first term represents the 
predicted signal for each data trace. Equation (1) assumes that the filter F s is exactly 
correct, so that this term is the exact signal component Xs/ of each data trace D,. 

The second term on the right hand side of Equation (1) is the results of the convolution 
of the swell noise with the signal filter, and this is termed the "residual". If the filter F s 
is exactly correct, the residual will not contain any of the signal of interest. 

In practice, the filter F s derived from the data traces D, is unlikely to be exact so that the 
residual of the right hand side of equation (1) after a single application of the filter will 
contain some components of the signal of interest It is therefore preferred to adopt an 
iterative approach, so that the residual term can be attenuated through further iterations 
thereby reducing the signal component of the residual. This is done by re-applying the 
filter F s to the full right-hand side of equation (1), and adding the results to the 
corresponding results obtained by the first iteration. The summation of the results of the 
two iterations will enhance the quality of the signal (the terms in Xs) and attenuate the 
noise (the term in X N ) through destructive interference. Further iterations may be 
applied as necessary, until the noise has been attenuated to a satisfactory level. 

Equation (1) is performed at each frequency in the swell noise range, using the 
appropriate F s (i.e. Fs(f L ) for data at frequency f L , Fs(f L +i) for data at frequency f L +i, 
etc). 
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If the lower frequency limit of the swell noise range, F L , is greater than zero, data at 
frequencies below F L will not be filtered. In this case, the output will consist of: 

a) for frequencies of F H +i or above - unchanged signal-only data (since there is no 
. swell noise in this range in the initial data); 

b) for frequencies in range from F L to F H - signal-only data (obtained by removal of 
swell noise); and 

c) for frequencies below F L - signal and other noise (not swell noise). 

If the lower frequency limit of the swell noise range, F L , is zero, the output will consist 
only of (a) and (b) above. 

Figure 6 is a schematic block diagram of an apparatus 11 that is able to perform a 
method according to the present invention. The apparatus 11 may be used in a marine 
seismic surveying arrangement in place of the conventional processor 8 of Figure 1, or 
it may be used to process pre-existing data. (If the apparatus 11 is used in a 
conventional marine seismic surveying arrangement of 1he type generally shown in 
Figure 1, the seismic surveying arrangement requires means for acquiring seismic data 
in the frequency-space domain and these are conveniently provided by receivers for 
acquiring seismic data in the time domain and transform means for transforming the 
data to the frequency-space domain. The transform means may be incorporated in the 
apparatus 1 1 or they may be additional to it.) 

The apparatus 1 1 comprises a programmable data processor 12 with a program memory 
13, for instance in the form of a read only memory (ROM), storing a program for 
controlling the data processor 12 to process seismic data by a method of Ihe invention. 
The apparatus further comprises non-volatile read/write memory 14 for storing, for 
example, any data which must be retained in the absence of a power supply. A 
"working" or "scratch pad" memory for the data processor is provided by a random 
access memory RAM 15. An input device 16 is provided, for instance for receiving 
user commands and data. One or more output devices 17 are provided, for instance, for 
displaying information relating to the progress and result of the processing. The output 
device(s) may be, for example, a printer, a visual display unit, or an output memory. 
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Sets of seismic data for processing may be supplied via me input device 16 or may 
optionally be provided by a machine-readable data store 18. 

The results of the processing may be output via the output device 17 or may be stored. 

The program for operating the system and for performing the method described 
hereinbefore is stored in the program memory 13, which may be embodied as a 
semiconductor memory, for instance of the well known ROM type. However, the 
program may well be stored in any other suitable storage medium, such as a magnetic 
data carrier 13a (such as a "floppy disk") or a CD-ROM 13b. 
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CLAIMS: 



1 . A method of processing seismic data, the data being in me frequency-space 

domain, the method comprising: 

detennining a prediction filter from the seismic data at a first frequency at which 

swell noise is not present; and 

applying the prediction filter to seismic data at a second frequency at which 
swell noise is present, thereby to attenuate swell noise in me seismic data at the second 
frequency. 

2. A method as claimed in claim 1 and further comprising the steps of: 
determining a prediction filter from the seismic data at a third frequency at 

which swell noise is not present, the third frequency not being equal to the first 
frequency; and 

applying the prediction filter to seismic data at a fourth frequency at which swell 
noise is present, the fourth frequency not being equal to the second frequency, thereby 
to attenuate swell noise in the seismic data at the fourth frequency. 

3. A method as claimed in claim 2 wherein the difference between the first 
frequency and the second frequency is equal to the difference between the third 
frequency and the fourth frequency. 

4. A method as claimed in claim 1, 2 or 3 and comprising, for each frequency in 
the range for which swell noise is present in the seismic data, determining an associated 
prediction filter from the seismic data at a frequency at which swell noise is not present. 

5. A method as claimed in claim 4 wherein the frequency range for which swell 
noise is present in the seismic data extends from a frequency f L to a frequency f H , and 
the associated prediction filter for frequency f s is determined from seismic data at a 
frequency fi-L+H+i- 
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6. A method of marine seismic surveying comprising actuating a source to emit 
acoustic energy; acquiring seismic data in the frequency-space domain; and processing 
the data according to a method as defined in any preceding claim. 

7. A method as claimed in claim 6 wherein the step of acquiring seismic data in the 
frequency-space domain comprises acquiring seismic data in the time domain and 
transforming the data to the frequency-space domain. 

8. A method as claimed in any preceding claim and comprising the further step of 
transforming the processed data to the time domain. 

9. An apparatus for processing seismic data, the data being in the frequency-space 
domain, the apparatus comprising: 

means for determining a prediction filter from the seismic data at a first 
frequency at which swell noise is not present; and 

means for applying the prediction filter to seismic data at a second frequency at 
which swell noise is present, thereby to attenuate swell noise in the seismic data at the 
second frequency. 

10. An apparatus as claimed in claim 9 and further comprising: 

means for determining a prediction filter from the seismic data at a third 
frequency at which swell noise is not present, the third frequency not being equal to the 
first frequency; and 

means for applying the prediction filter to seismic data at a fourth frequency at 
which swell noise is present, the fourth frequency not being equal to the second 
frequency, thereby to attenuate swell noise in the seismic data at the fourth frequency. 

11. An apparatus as claimed in claim 9 or 10 and adapted to determine, for each 
frequency in the range for which swell noise is present in the seismic data, an associated 
prediction filter from the seismic data at a frequency at which swell noise is not present 
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12. An apparatus as claimed in claim 11 wherein me frequency range for which 
swell noise is present in the seismic data extends from a frequency f L to a frequency fe. 
and the apparatus determines, in use, an associated prediction filter for frequency ft from 
seismic data at a frequency fi-L+m-i. 

13. An apparatus as claimed in any of claims 9 to 12 and comprising a 
programmable data processor. 

14. A seismic surveying arrangement comprising: a source of acoustic energy; 
means for acquiring seismic data in ihe frequency-space domain; and an apparatus as 
defined in any of claims 9 to 13. 

15. A seismic surveying arrangement as claimed in claim 14 wherein the means for 
acquiring seismic data in me frequency-space domain comprises one or more receivers 
for acquiring seismic data in the time domain and transform means for transforming the 
data to the frequency-space domain 

1 6. A storage medium comprising a program for a data processor of an apparatus as 
defined in claim 13. 

17. A storage medium containing a program for controlling a data processor to 
perform a method as defined in any of claim 1 to 5. 
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